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CN ■ Abstract 

We consider an equation 



Ph. 



, ^ - 2a 2/3 27 „ 

(U) = Uxx + Uyy + Uzz H Ux H Uy H U-^ — 

X y z 



in a domain Rg" = {{x,y,z) : x > Q, y > Q, z > 0}. Here a,/3,7 are constants, moreover 
< 2q:, 2/3, 27 < 1. Main result of this paper is a construction of eight fundamental solutions for 
above-given equation in an explicit form. They are expressed by Lauricella's hypergeometric functions 
of three variables. Using expansion of Lauricella's hypergeometric function by products of Gauss's 
hypergeometric functions, it is proved that the found solutions have a singularity of the order 1/r at 
. r ^ 0. Furthermore, some properties of these solutions, which will be used at solving boundary-value 

problems for afore-mentioned equation are shown. 
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\Q ' 1 Introduction 

' Known that fundamental solutions (FSs) have an essential role at studying partial differential equations 

(PDEs). Formulation and solving of many local and non-local boundary- value problems (BVPs) are based 
' on these solutions. Moreover, FSs appear as potentials, for instance, as a simple-layer and double-layer 

On . potentials in the theory of potentials. 

' Explicit form of FSs gives a possibility to study considered equation in detail. For example, in the 

works [3-5] by J.Barros-Neto and I.M.Gelfand, FSs for Tricomi operator, relative to an arbitrary point 
in the plane were explicitly calculated. Shown that found FSs clearly reflect the change of type of the 
Tricomi operator across the cc-axis. We also mention J.Leray's work [12], where a general method, based 
' upon the theory of analytic functions of several complex variables, for finding FSs for a class of hyperbolic 

linear differential operators with analytic coefficients was described. In particular, he showed how his 
method could be used to obtain, in the hyperbolic region, a FS for the Tricomi operator relative to a 
point (0,6). Among other results on this direction we would like to note work by M.Itagaki [11], where 
three-dimensional high-order fundamental solutions for modified Helmholtz equation were found. Found 
solutions can be applied with the boundary particle method to some 2D inhomogeneous problems, for 
example, see [8]. 

Various modifications of the equation 

iQ,/3,7 (m) = + Uyy + Uzz + — + — u,„ + —u^ = 0, < 2a, 2/3, 27 < 1, a, /?, 7 = const, (1) 

X y z 

in two-dimensional case were considered in many papers [1, 7, 13-18]. For example, in [7] the author is 
concerned with generalized bi-axially symmetric potentials (GBSP), i.e., with functions 

oo 

u(a;,y) = ^a„r2"P^"-« (cos 20), 

n=0 
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where x = rcos9, y = rsiuO, a,/3 > —1/2, Pn (cos 20) are Jacobi polynomials. These functions occur 
as solutions of the equation 

2/3+1 2a + 1 

X y 

Theorems on the convergence of above series and growth properties are obtained in terms of the coefficients 
In [1], the elliptic or ultra hyperbolic equation 




n s 

is considered, where the constants ai, «„, /3i, /Jg 7 are real parameters, = a;| ± ^ j/? and 

i=l i=l 

establishes some properties of the operator L. An example is the Kelvin principle for the equation 
(2): If u{x,y) is a solution of (2), then the function a = r~'^u{^,r]) is also a solution of the same 

n s 

equation, where = Xi/r^, rji = yi/r'^, (p = n + s — 2 + J2 + J2 Pi- For a class of iterated 

i=l i=l 

equations L^u = 0, it was shown that if Lu = 0, then [r^^u{x,y)\ = LP ['''^^"'^m 77)] = for 
j = 0, ...,p — 1; if, the function u {x,y) is a homogeneous solution of degree A of the equation Lu = 0, 
then LP [r^J-^'-^^^u {x, y)] = Lp [r20+^)u (^, rj)] = for j = 0, 1, ...,p- 1. 

FSs for elliptic equations with singular coefficients are directly connected with hypergeometric func- 
tions (HFs). Therefore, basic properties such as decomposition formulas, integral representations, formu- 
las of analytical continuation, formulas of differentiation for HFs are necessary for studying FSs. 

Since afore-mentioned properties of HFs of Gauss, AppcU, Kummer were well-known [2], results on 
investigations of elliptic equations with one or two singular coefficients were successful. Up to publication 
results on LauriccUa HFs by A.Hasanov and H. M.Sri vast ava [9, 10] there was no possibility to investigate 
FSs for the equation (1), since their consist LauriccUa HFs of three variables. 

In the present paper we construct eight fundamental solutions for the equation (1) in an explicit form. 
Using decomposition formulas for LauriccUa HFs of three variables by simple Gauss HFs, we prove that 
found FSs have a singiUarity of order l/r at r ^ 0. 

The plan of this paper is as follows. In Section 2 we briefly give some preliminary information, which 
will be used further. Also some constructive formulas for the operator iQ.fj,- are given. In section 3 we 
describe the method of finding FSs for considered equation and show what order of singularity will have 
found solutions. 



2 Preliminaries 

Below we give some formulas for Euler gamma-function. Gauss HF, Lauricella HF of three variables, 

which will be used in next section. 

Known that Euler gamma- function r(a) has properties ([6], pp. 17-19, (2), (10), (15)) 



T{a + m) = r(a)(a)„; T{a + 1/2) = r(l/2) = ^. 



(3) 



Here {a)m is a Pochgammer symbol, for which an cquahty {a)m+n = (a)?n(a + w)„ is true ([6], p. 67, (5)). 

^ (ci) (b^ 

A function 2^^1(0, b; c; x) — > ——^ — p-x" is known as Gauss HF and an equality 

^0 ('^)""' 

2Fi(a, b; c; 1) = ^l^^^^.^. ~ " ~ , c ^ 0, -1., -2, Re{c - a - 6) > (4) 
1 (c — 0)1 (c — a) 

is hold ([6[, p. 73, (14)). Moreover, the following autotransformer formula ([6], p. 76, (22)) 

2Fi(a, b; c; a;) = (1 - a;)"^ 2^^1(0 - a, b; c; -^—) (5) 

X — I 
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(6) 



is valid. 

The following system called as the system of hypergeometric equations of Lauricella ([2], p. 117): 

C (1 - - Cw«r, - CCt^cc + [2a - (a + /3 + 7 + i + a + l) - aijuj^ 

—a(uj(; — (a + /3 + 7 + i) ao; = 0, 

77 (1 - ry) - Cv^^^n - vCi^riC + [2/3 - (a + /? + 7 + | + /3 + l) ??] - P^^^i 

-/3Cw^ - (a + /3 + 7 + i) /3w = 0, 

C (1 - C)'^CC - ^Cwjc - ^Cwr,c + [27 - (a + /3 + 7 + i + 7 + 1) C] - 7^wj 

-jrioJri - (a + /? + 7 + 5) 70; = 0. 

The system (6) has the following solutions ([2], p. 118): 

wi(^,r?,C) = id'^ (|a + /3 + 7+^;a,/3,7;2a,2/3,27;^,r?,c) , (7) 

^^2 r?, C) = ) + /? + 7 + 1 1 - a, /3, 7; 2 - 2a, 2/3, 27; ^, 77, , 

W3 (C, ^, C) = V^-^f^F^^^ - /3 + 7 + I a, 1 - /3, 7; 2a, 2 - 2/3, 27; ^, 7?, , 
^^4 r?, C) = C'^^'J^r^ (^a + /3 - 7 + I a, /3, 1 - 7; 2a, 2/3, 2 - 27; ^, 7?, , (10) 
^5 r?, C) = e-^"v'-^^F^I^ (^-a - /3 + 7 + ^; 1 - a, 1 - /3, 7; 2 - 2a, 2 - 2/3, 27; ^, 7?, c) , (11) 
'^6(^,r?,C) =$'-'"C'-M'^ (^-a + /3-7+^;l-a,/3,l-7;2-2a,2/3,2-27;^,7?,c) , (12) 

(13) 



(8) 
(9) 



^7 (e, V, C) = 7?i-2/3ci-27F^3) _ /3 _ ^ + c,, 1 - /3, 1 - 7; 2a, 2 - 2/3, 2 - 27; ^, r,, c) , 

W8(^,»?,C) .^.^ 

^ ^i-2«^i-2/3^i-27j^^3) _ ^ _ ^ + 7. 1 _ 1 _ 1 _ ^. 2 _ 2a, 2 - 2/3, 2 - 27; r,, C) , ^ ^ 
where ([2], pp. 114 - 115 (1), (5)) 

r-,(3) / , , , \ \ ^ {Cl)i+j+k{bl)i{b2)j{b3)k i i k /I I 



p(3). , , . . r(ci)r(c2)r(c3) 

i^, (a; 6„ 6„ 63; c„ c„ 03; y, z) = r(,^)r(,^)r(63)r(c, - 6,)r(c. - 6.)r(c3 - 63) 
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yyy t?i-H^^-H^^-i(i-ti)^^-''^-i(i-t2)^='-'^-'(i-t3)^^-''^-'(i-xfi-yt2-2t3)-"dti«3, 



000 

iJe Cj > i?e 6i > 0, i = 1, 2, 3. 

The following decomposition formula ([9], p. 117,(14)) 



t;,(3) . u U h \ ST^ ('^)i+m+7i (^l);+m (^2);+„ (&3)to+„ 



z,„,„=o (c2)j+„ (c3)„+ J!m!n! ^^^^ 

X 2i^i (a + ; + 7n, 61 + ? + TTi; ci +? + Tu; a;) 2-Fi (a + i + m + n, 62 + ^ + «; C2 +/ + n; y) 
X2F1 (a + ? + 7n + 71, 63 + m + 7i; C3 + m + 7i; 2r) 



and formula of differentiation ([2], p. 19, (20)) 

-Ff^ {a + i + j + k,bi+iM + 3, bs + fc; ci + i, C2 + j, C3 + A;; x, y, z) 



d^+'+^Ff _ {a),^,+, (fei), (62), (63), ^(3) 



dx^yiz'' (ci) . {C2)j {03)^. 

(16) 
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are valid. 

One of the reasons why we consider the equation (1) is an existing of constructive formulas for the 

operator La.p.-y, which give a possibility to investigate the operator at various values of parameters a, /3, 7. 
Remark 1. The following constructive formulas 

ia,/3,7 (y'l^^uj = y^-_^^L^,i-p,^ {u) , (17) 

Lc,pn W-^'^z^-^'^u) = a;i-2"zi-27L^_^_^_^_^, , (is) 
ia,/3,7 {x'-'"y'-'''z^-'^u) ^ {u) (19) 

are true. 

Let us prove (19). Note, equalities (17), (18) can be proved similarly. First we find 
L^,0,, (x'-'-y'-'^z'-^^u) 

= (x'-^-y'-^^z'-^-'u)^^ + (xi-2"yi-2/3^i-27„)^^ + {x'-^'^y^-^^ z'-^-'u) 
X y ^ z ^ 

Since 

(3.1-2a 1-2/3^1-27^) 
V ^ /XX 

= ari-2«yi-2/3^i-27^^^ + 2(1 - 2a)x-2«yi-2/3^i-27„, - 2a(l - 2a)x-^"-^y^-^^ z^-^^u, 

(2.1-2a 1-2/3^1-27^) 

= a;i-2"yi-2/32i-27«^j, + 2(1 - 2/3)x^-2«y-2/32i-27i,^ - 2/3(1 - 2/3)xi-2"y-2/3-i2i-27 



(^l-2ayl-2/3^1-27y)^^ 

= + 2(1 - 27)a;^-2"yi-2/3^-27y^ _ 3^(1 _ 27)x^-2«yi-2/3^-27-iu^ 

^ = 2ax-2«yi-2/3^i-27„, + 2a(l - 2a)x-^''-^y^-^l^ z^-^^u, 

^ (a;i-2"yi-2/5^i-27„) = 2/3a;i-2«y-2/3^i-27y^ + 2/3(1 - 2(i)x^-^"y-^^-^z^-^''u, 

y ^ 

^ (xl-2"yl-2/3^1-27y) = 27xl-2"yl-2/5^-27^,^ + 27(1 - 27):cl-2"yl-2/5^-27-l^, 

after some simplifications we get 



^l-2ayl-2/3^1-27i^_^^^_^_^_^(^). 
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3 Finding fundamental solutions for the equation (1) 

Solutions of the equation (1) we shall search as follows 

u = P{r)uj{^,V,0, (20) 

where 



H = - xo) + (y - ya) + {z- zo) , 
rl = {x + xof + {y- Vof + {z- zof , 
r^ = {x- Xo) +{y + Vo) +{z- Zo) , 
rl = {x- xo)'^ + {y- yof + {z + zof , 



(21) 



^=^,.= ^,C=^, Pir)^irr"-'-''-'^. (22) 



Here V(a;, y, z) G Rg", but {xo, yo, zo) G is any fixed point. Substituting (20) into (1) we get 

+ A2U>riTi + ^S^CC + BiLJ^ri + -B2W|^ + ^3^,,^ + CiU)^ + C^Wj^ + + Duj = 0, (23) 

where 

Al=P [e +ey+ ] ,A2=P[VI+Vl+ V^] ,As=P [C + C'y + C] , 
Bl = 2P [^^r]^ + ^yfjy + ^,r],] , B2 = 2P ^^Cx + ^yCy + LCz] , S3 = 2P [r?^Cx + VyCy + VzCz] , 

C\ = P (^xx + ^yy + 6z) + 2 (Px^x + Py^y + Pz^z) + J' Kx — + ?y — + — 



C2 = P (r?xx + Vyy + Vzz) + 2 {PxVx + PyVy + PzVz) + P (vx— + Vy— + Vz—j , 

\ X y z J 

Cs = P (Cxx + Cyy + Czz) + 2 (PxCx + PyCy + PzCz) + P ( + Cy— + Cz — ) , 

\ x y z ) 

D = Pxx "H Pyy ~^ Pzz ~t~ Px "t~ Py ~t~ Pz • 

X y z 

After several evaluations we find 

= -4P (r^) x-^xo^ (1 - , ^2 = -4P (r2) y-^yoT) (1 - r?) , (24) 

A3 = -4P (r2) z-^zoC (1 - C) , Si = 4P (r2) x-^xoiv + 4S (r^) jy-^yoC??, (25) 

B2 = 4P (r^)"' x-^xoCC + 4.P (r^)"' ^-^zo^C, B3 = 4P (r^)"' y-^yoryC + 4P (r^)"' z-^zotjC, (26) 

Ci = -4P(r2)"\-ia;o [2q - (a + /3 + 7 + i + a + l) C] + 4P (r^)"' y-iyo/^C .o?! 
+4P(r2)-^-izo7?, ^ ^ 

C2 = -4P (r2) y-iyo [2/3 - (a + /3 + 7 + i + /3 + l) r?] + 4P (r^) x-ixoar? 
+4P(r2)"^z-izo7r?, 

C3 = -4P (r2)"' z-izo [27 - (a + /3 + 7 + i + 7 + 1) C] + 4P (r^) x-^aJoaC .39) 
+4P(r2)-\-iyo/3C, 

= 4P (r2) x-^xo (a + /3 + 7 + 5) a + 4P (r^) y-iyo (a + /? + 7 + ^) /? /ggx 

+4P (r2)"' z-izo (a + /3 + 7 + 5) 7- 

Substituting equalities (24) - (30) into (23) we obtain the system of hypergeometric equations of Lauricella 
(6), which has solutions (7) - (14). 
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Considering (7) - (14), from (20) we obtain eight FSs of the equation (1) 

91 {x, y, z; xq, yo, zo) = h (ry~^~^~'^ F^^^ (^a + p + -y + ^■,a, /3, 7; 2a, 2/3, 27; ^, rj, , (31) 

— / \ (32) 

= k2 (r2) 2 + ^ + ^ + _. 1 _ /j, ^; 2 - 2a, 2/3, 27; ^, r?, cj , 

93 (a;,y,^;;a;o,yo,2;o)^ 

= h (r2) 2 yi-^'^yi-^'^Fi^') _ ;3 + ^ + I c., 1 - /3, 7; 2a, 2 - 2/3, 27; V, c) , ^^^^ 

g4(a;,2/,2;;a;o,yo,^;o) 
3 

-a-/3+7-7r , o., / . . .3 _ „ \ (34) 



fc4 (r^) " 2 ^1-27^1-27^^3) + ^ _ ^ + 3 . 1 _ ^. 2c,, 2/3, 2 - 27; r?, c) , 



5 

q5{x,y,z;xo,yo,zo) = k5{r^) 2 



x^i-2ayi-2/3^i-2ayi-2/3^j3) _ ^ + ^ + 5 _ ^ _ ^ _ ^. 2 _ 2a, 2 - 2/3, 27; ^, r?, c) 



5 

q6{x,y,z;xo,yo,zo) = keir"^) 2 



X^l-2a^l-2^^1-2a^l-2T,^|3) + ^ _ ^ + ^ . ^ _ „^ j _ ^. 2 _ 2a, 2/3, 2 - 27; V, C) , 



5 

_Q|_|_^_|_.y_ — 

q7{x,y,z;xo,yo,zo) = kr (r^) 2 



xyi-2/5^i-27yi-2/3^i-27j,|3) _ ^ _ ^ + 5 . ^ _ j _ ^. 2a, 2 - 2/3, 2 - 27; ^, 7?, c) , 

7 

(a.,2/,^;a;o,yo,^o) = fcs (r^)"^^^"" 2 ^i-2aj^i-2/3^i-2T,^i-2a^i-2/3^i-2, 
xFi^) ('-a-/3-7+^;l-a,l-/3,l-7;2-2a,2-2/3,2-27;C,r?,C 



(35) 



(36) 



(37) 



(38) 



Here ki {i = 1,8) are constants. They will be determined when we solve boundary- value problems. 
Let us show that the found solutions (31) - (38) have a singularity. 

We choose a solution qi {x, y, z; xq, yo, zq). For this aim we use expansion for hypergeometric function 
of Lauricella (15). As a result solution (31) can be written as follows 

1 oc. (a + /3 + 7+J) (a),+„ (/3);+„ (7)„+„ 

qi{x,y,z;xo,yo,zo) = ki(r^) 2 > -— ^ ^ jr—^ 

i^^^o (2a);+„(2/3);+„(27)„+„Z!m!n! 

2 \ l-\-m / 2 \ / 2 \ "^^"^ 

2 X (39) 



X2Fi^a-|-/3-|-7-|-^-|-/-|-m-|-n,7-|-m-|-n;27-|-m-|-n;l - . 

Using formula (5), we rewrite (39) as 

gi (x,y,^;;a;o,yo,2;o) = fcir"^ (r?) " (rf) (r^) / (r^,rf ,ri,r|) , (40) 



/ 1 ■ J.2 

X2Fi a + /3 + 7+ -+ / + m,a + Z + m;2a + Z + m;l-4 
\ 2 

xsi^i (a + /3 + 7+^+Z + m + n,/3 + Z + n;2/3 + Z + n;l 



1 r2 
- + Z + m + n, /3 + Z + n; 2/3 + Z + n; 1 - -„ 

2 r-^ 

1 , ^ ^ T 
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where 



2 / 

f{r^,rl,rlrl)= 



2 / \ 

2 , , „ , 

/ 1 

X2-F1 a-/?-7--,a + / + m;2a + Z + m;l 5- 

V 2 rf 

X2F1 ^/3-a-7-^-m,/? + Z + n;2/3 + Z + n;l-^ 

/ 1 r 

X2-F1 l7-Q:-/J---Z,7 + m + n;27 + m + n;l ^ 

\ ^ ''3 

Below we show that / (r^, , r2, r^) will be constant at r ^ 0. 
For this aim we use an equality (4). Then we get 



1 

2' 

1' 



2-Fi { a — /3 — 7 — -, a + I + m;2a + I + m; 1 2 



r=0 



r(/3 + 7+2jr(2a + Z + m) 
r(a)r(a + /? + 7+ i+ Z + m 



Considering first formula in (3), we obtain 



1 

2F1 { a — (3 — 'y — -, a + I + m;2a + I + m;l j 

2 ^1 

r(/3 + 7+i)r(2a) (2a);+„ 



r=0 



r(a)r(a + /3 + 7+i) ( « + /J + 7 + 1 



l+m. 



Similarly we get 



1 

2F1 [ 13 - a - -y - - - m, (3 + I + n;2(3 + 1 + n;l - — 



r (2/3) r a + 7 



i) (2/3w(a + 7+i)^ 



r=0 



T{p)T[a + (3 + -i+- \ [a + (3 + ^+- 



2-fl 

'3 



r(27)r(^a + /3+i) (27W„ + ^ 



r=0 



r(7)r(a + 7 + /3+ij U + 7 + + 1 



Taking (42) - (44) into account from (41) at r ^ we have 



r(2a)r(2/3)r(27)r(a + /3+i )r(a + 7+^ )r(/3 + 7 + ^ 



r(a)r(/3)r(7)r3(a + /3 + 7 + l 



E 



a + 7 + (a + /? + ^ (a),+„ (/3);+„ (7)„+„ 



(45) 



!,m,n=o / a + 7 + /3 + i J (a + /3 + 7+^) Z!m!n! 

We use formulas in (3) to simplify (45). Since 



E 



l-\-m 



l-\-m 

I 



n=o(a + 7 + /J+ij (a + /3 + 7+^) l\m\n\ 



E 



/ /+m+n V / l-\-m 

(^a + 7 + ^) (a + /3 + ^ {a\ (a + l)^ (/3), {(3 + 0„ (t)^ (t + m\ 



!,m,n=0 ,^.1,,. \ I f . ^ . .1\ 



2 



Q; + 7 + /3+-+Z + m a + /? + 7+ - l\m\n\ 



2 /\V 2/ 



l-\-m J 



^ V + g (/^ + 0„(7 + m)„ 



^^a + /3 + 7+-') 1 ;,™=oU + 7 + /J+l+/ + m) n! 



2 

1 



r a + /3 + 7+- r a , , ^ / a ^ ^ ^ , n 

V 2/ V ) ^ Wlif^)l {Vrn(a + l)r 

r ( a + 7 + i ) r ( a + /? + i ) i=o ( a + /3 + 7 + i ) U m=o I a + P + ^ + - + I ] m\ 



r2(a + /3+ijr(^a + 7+ i)r(/3 + 7+ ^)i=o(a + /3 + 7 + 
r2 ( a + 7 + /3 + i 



r(a + /3+i)r(a + 7 + ^)r('^ + 7 + ^ 



2 

from (45) follows 



r (a) r (/3) r (7) r a + + 7 + 



Expressions (46) and (40) gives us possibility to conclude that the solution gi (x, y, 0; xo, yo, ^0) re- 
duce to infinity of the order at r — > 0. Similarly it is possible to be convinced that solutions 

qi {x, y, z: xn. yo, zq) , i = 2, 3, 8 also reduce to infinity of the order when r — > 0. 
It is directly checkable that constructed functions (31) - (38) possess properties 

ar^^AgJ 0, y^^^qi\ , = 0, z^-^^qA , = 0, 
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Let be convinced that 



d 

a;'"^57L=o = 0' '^4=0 = 0' «4=o = 0' 

^'"#:'?i|.=n = 0- (47) 



Using formula of differentiation for Lauricella HF of three variables (16) we have 



X 



91 {x, y, z; xo, yo, zq) = -2{x - xq) (^a + (3 + ^ + ^ h (r^) 



2\ -a-/3-7- 



dx 

(3) 



F^i^^ ( a + /? + 7 + I a, /3, 7; 2a, 2/3, 27; ^, r?, C 



^F;^-^^ (a + /3 + 7 + -;a + l,/3,7;2a + l,2/3,27;C,r7,C 



+ (a + /3 + 7+2;"'/3 + l>7;2a,2/3+l,27;^,r/,C 



^^f) (a + /3 + 7+|a,/3,7 + l;2a, 2/3,27+ l;C,r?,C 



-4xo(^a + /3 + 7+^jfci(r2) " ^ ^ ^ (^a + /3 + 7 + ^; a,/3,7; 2a, 2/3, 27; ^, 7?, C 

From th{^ last equality one can easily be convinced that (47) is true. 

Remark. Afore-mentioned properties of found fundamental solutions will be used at solving various 
boundary-value problems for the equation (1). For instance, fundamental solution (31) with its properties 
will he used at solving problem N for the equation (1). 
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